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Abstract 

An analysis of the semiclassical regime of the quantum-classical transition is given for open, 
bounded, one dimensional chaotic dynamical systems. Environmental fluctuations - characteristic 
of all realistic dynamical systems - suppress the development of fine structure in classical phase 
space and damp nonlocal contributions to the semiclassical Wigner function which would otherwise 
invalidate the approximation. This dual regularization of the singular nature of the semiclassical 
limit is demonstrated by a numerical investigation of the chaotic Duffing oscillator. 
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The process by which a classical dynamical system emerges as a sufficient 
approximation to a quantum dynamical system has been a major topic of dis- 
cussion since the inception of quantum mechanics. The singular nature of the 
semiclassical limit (h 0) [Ij lies at the center of this debate. Classical be- 
havior cannot emerge as the smooth limit of a closed quantum system with a 
nonlinear Hamiltonian as these classical evolutions violate the unitary symme- 
try of quantum mechanics [2]. Moreover, the symplectic geometry of a chaotic 
classical phase space generates infinitely fine structures which the uncertainty 
principle prevents a quantum dynamical system from tracking as t — oo 
Thus the pathologies associated with the semiclassical limit are most dramatic 
in classically chaotic dynamical systems. In these systems, small h expansions 
of physical averages have been shown to fail at a finite time This work ad- 
dresses how the incompatibility of quantum dynamics and classical chaos can be 
resolved in the framework of open quantum and classical systems. In particular, 
the stability of the semiclassical limit is recovered by an environment-induced 
coarse graining of both classical and quantum dynamics, from which a threshold 
condition is derived. 

As all realistic systems interact with their environment, the modern approach to un- 
derstanding the quantum-classical transition (QCT) relies on the open system paradigm. 
In this paradigm, the dynamical system is not considered in isolation but analyzed taking 
its external interactions into account. These environmental interactions are of two types 
depending on whether it is possible to make measurements on the environment or not. If 
the environment is in principle unobservable, then the system is described by the reduced 
density matrix obtained from the full system-environment density matrix and tracing over 
the environment. If, on the other hand, certain measurements are possible on the environ- 
ment, then the resulting reduced density matrix for the system depends on the results of 
these measurements. System evolution in this second case is therefore said to be conditioned 
o. tKe obse^fo. .e.^.ts g Q. The condit.o.ed syste,n state, as ,t evolves, is said to 
define a quantum trajectory and inequalities governing the existence of a classical trajectory 
limit of quantum trajectories in continuously measured quantum systems have now been 
obtained 0. Since it yields effectively classical trajectories - from a single realization of 
a measurement record - we will call this pathway the strong form of the quantum-classical 
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transition. 



If the environment is not amenable to observation, or if one decides to throw away the 
results of measurements on the environment - which amounts to the same thing - then the 
evolution of the reduced density matrix of the system is given by an unconditioned master 
equation, as one must average over all possible measurement outcomes. It is now no longer 
possible to obtain the classical trajectory limit as discussed above. One must now compare 
quantum and classical distributions (or, equivalently, the underlying moment hierarchy) 
against each other: This constitutes the weak form of the quantum-classical transition. For 
any given situation, if the inequalities defining the strong form of the QCT are satisfied, then 
the weak form of the QCT immediately follows. The reverse is not true, however. Finally, 
it should be noted that there is no direct connection between the inequalities defining the 
strong form of the QCT and the ones defining the weak form derived here. 

The weak version is just another way to state the conventional decoherence idea j^; how- 
ever, as discussed in Ref. mere suppression of quantum interference does not guarantee 
the QCT even in the weak form. To address this problem, our purpose here is to present a 
semiclassical analysis of the QCT for bounded, classically chaotic open systems focusing on 
the regularization of the singular — > limit via the weak form of the environmental inter- 
action. This is distinct from the state localization characteristic of the strong form of the 
QCT. Given a small, but finite, value of we aim to establish the existence of a timescale 
beyond which the dynamics of open quantum and classical systems becomes statistically 
equivalent if the environmental interaction is sufficiently strong. 

It is important to keep in mind that the results obtained in this work are intrinsically 
different from the inequalities describing the strong form of the QCT given in Ref. Ql- The 
analysis here refers to the coarse-grained distribution function (averaging over all measure- 
ments), whereas the analysis in Ref. Ql refers to the fine-grained distribution for a single 
measurement realization. The strong form of the QCT requires that the quantum distribu- 
tion be localized, and this is enforced by the conditions, for weak and strong nonlinearity, 
respectively Q]: 



In addition, the key low-noise condition on the trajectories is enforced by a double-sided 
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inequality _ 
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In the equations above, where k is the measurement strength, r] is the measurement 
efficiency, m is the particle mass, F is the force on the system measured at the expectation 
value (X) = X and s is the typical action of the system in dimensionless units. This second 
inequality first states that if the measurement is too weak, the measurement signal is noisy as 
it results from sampling a wide distribution. As the measurement strength is increased, the 
distribution becomes sharply peaked and there is very little noise in the measurement result. 
However, if the measurement is too strong, the resulting backaction drives the dynamics too 
hard and the trajectory becomes noisy. Either of these limits need not prevent a weak 
QCT from existing as the weak version of the QCT is not a direct consequence of state 
localization and therefore can exist even when these inequalities are not satisfied. It does 
not matter if the coarse-grained distribution is too wide, as long as the classical and quantum 
distributions agree, and, even if the backaction noise is large, the coarse-grained distribution 
remains smooth and the weak quantum-classical correspondence can still exist. 

We demonstrate that, for a bounded open system with a classically chaotic Hamiltonian, 
the weak form of the QCT is achieved by two parallel processes. First, the semiclassical 
approximation for quantum dynamics, which breaks down for classically chaotic systems 
due to overwhelming nonlocal interference, is recovered as the environmental interaction 
filters these effects. Second, the environmental noise restricts the foliation of the unstable 
manifold, the set of points which approach a hyperbolic point in reverse time, allowing the 
semiclassical wavefunction to track this modified classical geometry. Our approach explicitly 
incorporates both the stretching and folding typical of hyperbolic regions and the role of the 
environment as a filter on a phase-space quantum distribution. Thus our results are different 
from the purely local analysis of Ref. 0] and the study of decay of off-diagonal coherences 
in Ref. which focuses on the filtering aspect of the QCT. 

We examine a simple model of a quantum system weakly coupled to the environment so 
as to maintain complete positivity for the subsystem density matrix, p(t), while subjecting 
it to a unitarity breaking interaction. These conditions mathematically constrain the master 



equation to be of the so-called Lindblad form 



If this environmental interaction couples 



4 



to the position, as is often the case, the master equation will take the form: 

^ = -^[if,/5]-A;e.4^,[^,P]], (3) 

'or environmental coupling strength k^nv, subsystem Hamiltonian H and position operator X 
^, Q, Q| . We have neglected the dissipative environmental channel and kept the diffusive 
channel for two reasons: (i) the coupling to the environment is always assumed to be weak 
and the dissipative timescales are hence very long, longer than the dynamical timescales 
of interest, (ii) the weak form of the QCT arises only from the diffusive channel, hence 
dissipative effects are not of interest here. 

The master equation is often examined by taking the Wigner transform of the density 
matrix in its position representation. This yields the Wigner function, which is a quasiprob- 
ability distribution over phase space whose evolution can then be compared to the corre- 
sponding classical phase space distribution function Q]. The Wigner representation of the 
density matrix is given by 

Uq.V.t) = — - / rfXe-^^/Vg + (4) 
Znii J~oo L I 

In this representation, the master equation for the open Wigner evolution becomes 

= Lafn, + Lgf^ + D-^, (5) 

where the classical Liouville operator Ld = —{p/m)dx + {dxV)dp and the quantum Li- 
ouviUe operator = E„>i(^^"(-l)"/(22"(2n + 1)!)) 9f and D = Ti'^Kuv When 
Lg = 0, this equation reverts to the classical Fokker-Planck equation. It is important to 
keep in mind that the specific form of the diffusion coefficient (alternatively, kf,nv) depends 
strongly on the physical situation envisaged. Thus, if the master equation describes a we akly 
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coupled, high temperature environment, D = 2m'ykBT (7 is the damping coefficient) 
whereas for a weak, continuous measurement of position, the diffusion due to quantum 
backaction is D = %^k P]. The results in this paper hold for all of these cases. 

The closed-system evolution (fce™ = 0) of the quantum and classical Liouville {Lq = 0) 
equations for nonlinear systems yield very different results. The Wigner evolution is quickly 
dominated by nonlocal interference effects and has rapidly oscillating positive and negative 
components, whereas the classical evolution is always positive but develops very fine-scale 
structure. Numerical studies for a class of classically chaotic Hamiltonians carried out for 



small, but finite kenv, show that diffusion rapidly filters fast, nonlocal interference terms. 



15l llG^ . Thus, the QCT connects the 



but it also filters the fine-scale classical structure 
quantum master equation to the classical Fokker-Planck equation and not to the classical 
Liouville equation. Establishing the general conditions under which this occurs is the main 
aim of this paper. 

Once the QCT occurs, the effects of Lq in the evolution specified by Eqn. (0) are sub- 
dominant. Therefore, to understand how environmental noise acts in this limit, it suffices 
to consider the behavior of the corresponding classical Fokker-Planck equation. To do this, 
it is convenient to examine the underlying Langevin equations for noisy trajectories that 
unravel the evolution of the classical distribution function when Lg = 0. These are given by: 

dq = —dt (6) 
m 

and 

dp = f{q)dt + V2DdW (7) 

where f{q) = —dV{q)/dq and dW is the Wiener measure [{dWf = dt]. For constant D, 
one can just as well write \/2DdW = ^{t)dt, where ^{t) is a rapidly fluctuating Langevin 
force satisfying (^(t)) = and {C,{t)^{t')) = 2D6{t — t') over noise averages. To examine the 
effects of this force on the formation of the stable and unstable manifolds associated with a 
chaotic system we perform a perturbative expansion about a hyperbolic fixed point [q^q, 0), 
where f{qeq) = 0, using small- noise perturbation theory (i.e., assuming the effects of noise 
forces are small relative to the systematic forces on dynamically relevant timescales) jl7|. 
In the leading order approximation, we can separate the dominant systematic components 
from the noisy components via q(t) ~ qc(t) + QNit) and p(t) ~ Pc(t) + PN(t), leading to 
the usual Hamilton's equations for qc and pc, and to the coupled equations dq^^ = p^dt/m 
and dpN = mX^q^dt + ^dt, where mA^ = df{qeq)/dq defines the local Lyapunov exponent, 
A. These have the solution jisj ]. 

q{t) = qe, + C+e^* + C^e"^* 

with an analogous expression for p{t). Position and momentum are then rescaled according 
to q' = V Xmq and p' = p/y/Xm. This gives both dynamical variables the identical units of 
the square root of phase space area and renders the stable and unstable directions orthogonal. 



Projecting the solutions for g' and p' along the stable (-) and unstable (+) directions, we 
find, 

= V2A^C±e±^* ± -jL= f rf<(M)e±^(*-"). (9) 
V2Am "'o 

The effects of noise can now be analyzed on the distribution functions generated by the noisy 
trajectories. The average over all noisy realizations of the displacement in each direction 



is given by {u±) = \J IXmC ±e^^^ ^ as expected from a perturbation in the neighborhood 
of a hyperbolic fixed point. When examining the second order cumulants ones sees that, 
while the stable and unstable directions have variances of ±2^^(e^^^* ~ 1); the off-diagonal 
cumulant is — {u^){u^) = —Dt/{mX). In forward time, where the evolution of a 

trajectory is determined by the unfolding of the unstable manifold, this indicates that, as 
the trajectory evolves, it will simultaneously smooth over a transverse width in phase space 



of size \jDt/{m\). 

The above analysis implies a termination in the development of new phase space struc- 
tures at some finite time t*, whose scaling behavior can be determined, which need not be 
true in a non-compact phase space. The average motion of a trajectory is identical to its 
deterministic motion, so at time t, if the initial length in phase space is mq, its current length 
will be approximately mqc'^* as its forward time evolution will be dominated by its component 
in the unstable direction. Here A is the time-averaged positive Lyapunov exponent. If the 
region is bounded within a phase space area A, the typical distance between neighboring 
folds of the trajectory is given by 

l{t) ^ -e-~^\ (10) 

Mo 

where one should recall that lit) still carries the units of the square root of phase space area. 
However, since phase structures can only be known to within the width specified above, the 
time at which any new structure will be smoothed over is defined by 



V mX 

The above two equations can be used to determine t*, which is only weakly dependent on 
D and the prefactor in Eqn. ()10|) . Due to the smoothing, one does not see an ergodic phase 
space region, but one in which the large, short-time features which develop prior to t* are 
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pronounced and the small, long-time features which develop later are smoothed over by the 
averaging process. Therefore, to approximate noisy classical dynamics, a quantum system 
need not track all of the fine scale structures, but only the larger features which develop 
before the production of small scale structures terminates. (This classical suppression of 
structure parallels recent studies of chaotic advection-diffusion, where the efficiency of mix- 
ing is suppressed There, as here, the microscopic dynamics generates structure on 
increasingly smaller scales until diffusion terminates further development.) 

The conditions under which quantum dynamics can track this modified phase space ge- 
ometry can be examined by looking at the semiclassical Wigner function in the presence of 
an environmental interaction. As has been appreciated for a long time, the Wigner represen- 
tation provides a useful tool for semiclassical analysis 2^. By examining the semiclassical 
Wigner function, rather than the semiclassical wavefunction, the breakdown of the semiclas- 
sical approximation for chaotic systems can be clearly associated with a geometric picture. A 
general mixed state is an incoherent superposition of pure state Wigner functions, where an 
individual semiclassical pure state Wigner function can be formed by substituting the usual 
Van-Vleck semiclassical wavefunction in Eqn. (^. If we allow q to be perturbed by noise 
we can rewrite the classical action S{q,t) ^ S{qc,t) — \p2D Jq dt^{t)qc{t), as in Ref. jl^ . 
Following Berry we rewrite the action for the i-th solution to the Hamilton- Jacobi 
equation as S^iqc,t) = /^^^(g dq'pi{q',t) - dt'iJ(gc(0),p.(gc(0), t') {= dt'Hiit')) , where 
Pi{q,t) is the i-th branch of the momentum curve for a given q. If we average over all 
noisy realizations, after separating the contributions from identical branches, the following 
suggestive expression for the noise averaged semiclassical Wigner function obtains: 

1 /-"^ / DtX^^ 



271% J-oo 



(iXexp 



2r 



exp 



h 



/ dq'pi{q',t)-pX 
J a- 



+ 



2i Jij sin 

i<j 
ft 

dt (Tii — Ti. 



dq'pi{q',t) - / dq'pj{q',t) 
9c(0) JqciO) 



Qiq+,t)C,iq^,t) 



^\Mq+,t)\\J,iq.,t)\ 

for Jacobian determinant Ji{q,t) and transport coefficient Ci{q,t); q± = g± y and 



(12) 
(13) 



where is the i-th Maslov index 
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One can analyze the dominant contributions to the above integrals using the stationary 
phase approximation P]. If D = 0, these would contribute phase coherences at values of X 
which satisfy Pi{q + X/2, t) — X/2, t)—2pX = for the first term in the summation and 

Pi{q + X/2, t) + pj{q — X/2, t) — 2pX = for the second term, the former being the famous 
Berry midpoint rule. For a chaotic system, Berry argued that, due to the proliferation 
of momentum branches, Pi{q,t), arising from the infinite number of foldings of a bounded 
chaotic curve as t —>■ oo, a semi classical approximation would eventually fail since the 
interference fringes stemming from a given pi could not be distinguished after a certain time 
from those emanating from the many neighboring branches j^. While the precise value of this 
time has since been challenged numerically, the essential nature of this physical argument 
has remained valid js^. However, the presence of noise acts as a dynamical Gaussian 
filter, damping contributions for any solutions to the above equation which are greater than 
X ^ h/^/Dt. In other words, noise dynamically filters the long "De Broglie" wavelength 
contributions to the semiclassical integral, the very sort of contributions which generally 
invalidate such an approximation. If we rescale the above result and combine it with our 
understanding of how noise effects classical phase space structures, we can qualitatively 
estimate whether or not a semiclassical picture is a valid approximation to the dynamics. 
As discussed earlier, t* is the time when the formation of new classical structures ceases and 
l{t*) is the associated scale over which classical structures are averaged. The key requirement 
is then that the semiclassical phase filters contributions of size 

^ < lin (14) 



In other words, for a given branch, the phases with associated wavelengths long enough to 
interfere with contributions from neighboring branches will be strongly damped, and the 
intuitive semiclassical picture of classical phase distributions decorated by local interference 
fringes will be recovered. 

The weak form of the QCT is completed when the inequality (fT^ is satisfied. Substituting 
the scale of classical smoothing (fTT|) in this inequality, we find 

Dt* > Xmh. (15) 

[Note that the purely classical quantity t* is first independently determined by solving 
Eqn. (fTT|) and then compared to the RHS of the above equation.] While, the lefthand side 



of the inequality contains the mutually dependent t* and D, the right hand side depends 
only on fixed properties of the system and h. This condition therefore defines a threshold 
at which the semiclassical approximation becomes stable and which may be set in terms of 
either D ot t*. Once the threshold is met, t* becomes the time beyond which the semiclas- 
sical description is valid. The semiclassical nature of this condition becomes more evident 
on defining S = l{t*Y which, given that is an areal scale in phase space for the diffusion 
averaged dynamics, has dimensions of action. A physical interpretation is more apparent on 
rewriting (fT3j) as S = l{t*Y > which is readily identified as the usual condition for the 
validity of a semiclassical analysis. 




FIG. 1: Sectional cuts of Wigner functions (dashed lines) and classical distributions (solid lines) 
for a driven Duffing oscillator, after 149 drive periods, taken at p = for (a) D = 10~^; (b) 
D = 10^^; (c) D = 10^^. Parameter values are as stated in the text; the height is specified in 
scaled units. 

To illustrate these mechanisms we considered the Duffing oscillator with unit mass: H = 
/ 2 + Bx^ — Ax^ + Kx cos{ujt). This system was studied for a set of parameters [A = A = 10, 
B = 0.5 and u = 6.07) where the system is strongly chaotic (A = 0.57) and the 

dynamical evolution of its bounded motion is dominated by the homoclinic tangle of a 
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single hyperbolic fixed point. As a result, the long-time chaotic evolution can be completely 
characterized by the unstable manifold associated with that fixed point 2^. We chose 
= 0.1 in these calculations. 

The evolution of the Duffing system was calculated for both the classical and quan- 
tum master equations. Figure [T] shows sectional cuts at p = of the quantum and 
classical phase distribution functions for three different values of the diffusion coefficient, 
D = 10~^, 10~^, 10~^, after time T = 149 evolution periods. As t* varies slowly with D 
[Eqns. (jlOHlH) ]. in the three cases shown, t* ranges only from ~ 20 — 14 (note that t* <^T). 
It is easy to check that the inequality (fTH|l is strongly violated for D = 10~^, mildly vio- 
lated for D = 10^^, and approximately satisfied for D = 10^^. For D = 10^^, the classical 
and quantum sections show no similarities, as expected. The quantum Wigner function 
also shows large negative regions, reflecting strong quantum interference. On increasing D 
to 10"3 the magnitude of quantum coherence decreases dramatically and the classical and 
quantum slices have the same average value, as well as specific agreement on some large scale 
features. The two disagree, as expected, on the small scale structures. This indicates that, 
while the quantum and classical distributions do not exactly match, the Wigner function 
has now become sensitive to the larger features of the noise averaged classical distribution 
function, indicative of the transition to a semiclassical regime. At D = 10-2, there IS near 
perfect agreement between classical and quantum distribution functions, save on the small- 
est scales. When D is of order unity, the inequalities enforcing the QCT at the level of 
individual trajectories Ql are satisfied and the agreement is essentially exact. However, as 
indicated by Fig. HI^c), detailed agreement for quantum and classical distribution functions 
can begin at much smaller values of the diffusion constant. 

For more detailed evidence that, at D = 10~^, one is entering a semiclassical regime, 
in Fig. 121 we superimpose an image of the large scale features of the classical unstable 
manifold on top of the full quantum Wigner distribution a.t D = 10^^ after 149 drive 
periods [case (b) of Fig. 1]. The quantum phase space clearly exhibits local interference 
fringes around the large lobe-like structures associated with the short-time evolution of the 
unstable manifold. The appearance of local fringing about classical structures is direct 
evidence of a semiclassical evolution, where interference effects appear locally around the 
backbone of a classical evolution. This is in sharp contrast to the global diffraction pattern 
seen for D = 0, where the contributions from individual curves cannot be distinguished. 
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FIG. 2: Phase space rendering of the Wigner function at time t = 149 periods of driving. The 
early time part of the unstable manifold associated with the noise-free dynamics is shown in blue. 
The value oi D = 10^^ is not sufficient to wipe out all the quantum interference which, as expected, 
is most prominent near turns in the manifold. 

suppressing the appearance of any classical structure [l^. 

The arguments presented here do not apply directly to systems where the manifolds are 
non-compact (analogous to open flows in fluids). In this case noise does not necessarily 
lead to a termination of the structure production in phase space, hence a semiclassical 
approximation eventually fails. Thus the topological structure of the classical manifolds 
likely provides the physical underpinning of a proposed classification scheme for dynamical 
systems based on quantum-classical correspondence, originally stated in terms of spectral 
analysis of the density matrix jj]. An elaboration of both the mathematical analysis and 
numerical results will be presented in a companion paper {2^. 

Numerical simulations were performed on the Cray T3E and IBM SP3 at NERSC, LBNL. 
B.D.G. and S.H. acknowledge support from the LDRD program at LANL. The work of B.S. 
and B.D.G (partially) was supported by the National Science Foundation grant ^^0099431. 
This research is supported by the Department of Energy, under contract W-7405-ENG-36. 
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